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Abstract 



We demonstrate with a minimal example that in Filippov systems (dynamical systems 
governed by discontinuous but piecewise smooth vector fields) stable periodic motion 
with sliding is not robust with respect to stable singular perturbations. We consider 
a simple dynamical system that we assume to be a quasi-static approximation 
?H of a higher-dimensional system containing a fast stable subsystem. We tune a 

system parameter such that a stable periodic orbit of the simple system touches the 
i_i discontinuity surface: this is the so-called grazing-sliding bifurcation. The periodic 

^ orbit remains stable, and its local return map becomes piecewise linear. However, 

^ when we take into account the fast dynamics the local return map of the periodic 

orbit changes qualitatively, giving rise to, for example, period-adding cascades or 
^ small-scale chaos. 

^ Key words: vector fields with sliding, singular perturbation, discontinuity induced 

^ bifurcation 
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1 Introduction 



Filippov systems are dynamical systems governed by discontinuous but piece- 
wise smooth ordinary differential equations (ODEs) and they occur, for example, 
in the modelling of mechanical systems (with dry friction) [14], of electrical 
systems with switching [1], or of population dynamics (when, depending on 
the population size, either the individuals switch habitat or diet, or harvesting 
becomes restricted [3]). The phase space of a Filippov system is partitioned 
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Fig. 1. A trajectory x{t) through a point xq on the switching manifold T-Ls generically 
either (a) crosses the switching manifold T-Lg if the vector fields /+ and /_ both point 
in the same direction relative to T-Lg-, or (b) slides along T-Lg if /+ and /_ both point 
toward Hg- 

into domains where for each domain a different ODE governs the dynamics. In 
the simplest case we have two domains such that 



X 



iih{x)<Q, 
f+{x) if h{x) > 0. 



(1) 



The boundary between the domains is called the switching manifold: T-Lg — 
{h{x) = 0}. 

A special feature of Filippov systems is the so-called sliding mode, which means 
that a trajectory of the system does not follow any of the ODEs governing the 
domains but it rather 'slides' along the switching manifold, following a convex 
combination of the ODEs governing the adjacent domains: 



X 



fs{x) 



[dh{x)f4x)] ■ Ujx) - [dhix)Uix)] ■ U{x) 
dh{x)[U{x) - U{x)] 



(2) 



where '9' denotes partial differentiation (in the current case with respect to the 
variable x). Sliding occurs in all points Xq of the switching manifold Tig where 
both vector fields point toward Hg, that is, dh{x)f+{x) < and dh{x)f^{x) > 0. 
Figure 1 illustrates the two typically observed cases of interaction between 
flows and switching manifold. 

An important question for modelling is how sliding in a Filippov system is 
aff'ected by perturbations. If we add a small perturbation to /_ or or to the 
switching decision h (the derivative of the perturbation is also assumed to be 
small) then any exponentially stable periodic orbit or equilibrium of a Filippov 
system persists [4] and remains stable. This also applies to pseudo-equilibria 
(equilibria of the sliding flow fg, sitting exactly on the switching manifold) 
and to periodic orbits that have sliding segments. This persistence mirrors 
the results of classical bifurcation and invariant manifold theory for smooth 
dynamical systems [6]. 



2 



Another typical perturbation arising in the modelling process are stable singular 
perturbations. In a simple model one has replaced rapidly converging parts 
of the dynamics with their equilibrium, making the assumption that this 
equilibrium follows the slow dynamics quasi- statically. In a more complex of the 
same system (or in reality) the equilibrium of the fast dynamics is not attained 
perfectly, which constitutes a small perturbation. Practical examples of this type 
of perturbation are small capacitances and inductances in electrical circuits, 
imperfect rigidity in mechanical systems, or fast chemical reactions (or other 
processes) in biological systems. Again, for smooth dynamical systems classical 
theory [6] proves that all hypberbolic equilibria, periodic orbits and, more 
generally, normally hyperbolic invariant manifolds persist. That is, for example, 
an exponentially stable equilibrium or periodic orbit (and any of its bifurcations) 
observed in a simple model obtained by making quasi-static assumptions is also 
present when the fast dynamics is taken into account as long as the difference 
in time scale is sufficiently large. In general, in smooth dynamical systems any 
phenomenon that persists under regular perturbations (perturbations of the 
right-hand-side) also persists under stable singular perturbations. Fenichel's 
Theorem reduces hyperbolic singular perturbations to regular ones by proving 
the existence of a normally hyperbolic invariant manifold [6] . 

In order to find general statements how the dynamics of Filippov systems is 
affected by stable singular perturbations one has to study slow-fast systems of 
the form 

f/-(^'?/'^) if /i(a;, I/, £) < 0, 
\f+{x,y,e) a h{x,y,e) >0, 
£y= 9(x,y,e). (4) 

In system (3), (4) e measures the difference in the time scales between the 
evolution of the slow variable x and the evolution of the fast variable y. We 
assume that for £ = the (then algebraic) equation (4) can be solved for y for 
all X, and that this solution yo{x) satisfies the stability condition 

> — c > Respecd2g{x,yo{x),0) (5) 

with a uniform decay rate c (subscript 2 in the differential operator indicates 
differentiation with respect to the second argument). Condition (5) means that 
yo{x) is a locally exponentially stable equilibrium of the fast subsystem (4) if we 
treat the variable x in (4) as a parameter and set £ to in the right-hand-side 
of (4). 

Let us assume that the quasi-static approximation (called the reduced system 
from now) 

^ ^ if-{x, yo{x), 0) ff h{x, yo{x), 0) < 0, 

1 f+{x, yoi^), 0) ff h{x, yoix), 0) > 
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has a stable periodic orbit x{t) which is switching in xq = x{tQ) from the subdo- 
main {x : h{x,yQ{x),0) > 0} to sliding inside the manifold {x : h{x,yQ{x),0) = 
0} (as shown in Figure 1(b)). What happens to this periodic orbit if we include 
the singular perturbation effects by changing £ to a positive value? 

The references [9,10] studied this question under the simplifying assumption 
that the switching decision does not depend on the fast variable y (that is, 
= in (3), (4)) and proved that stable periodic orbits persist when e 
becomes positive but they acquire small boundary layers after switching. In 
practice, it is often difficult to check or guarantee the condition = 0. 
For example, a switch in an electronic circuit may depend on the voltage 
potential at a node j which is only a function of neighboring node potentials, 
which (let us assume) are all slow variables. However, taking into account a 
small parasitic capacitance (of order e) affecting node j will change its voltage 
potential into a fast variable, thus, making the switching decision dependent on 
a fast variable. The same effect occurs more generally in modelling: one often 
studies the simple model {e — 0) only, and it is impossible to tell if switching 
decisions depend on fast variables for any of the possibly significant singular 
perturbations without studying the more complex model (e > 0). This means 
that the assumption = of [9,10] is mathematically convenient because it 
allows to prove persistence of stable periodic orbits but it does not cover all 
cases of practical interest. 

Our paper studies the problem how the dynamics near a periodic orbit changes 
under a stable singular perturbation if the switching decision depends also on 
fast variables: ^2/1 7^ in (3). It turns out that there are qualitative changes if 
we increase e from to a positive value. Possible scenarios are, for example, a 
period-adding cascade on a parameter range of order e, or small-scale chaos 
(also of order e) around a periodic orbit that is exponentially stable for £ = 
but is uniformly unstable for all £ > 0. 

The paper is outlined as follows. In Section 2 we construct a simple example that 
shows how arbitrarily small stable singular perturbations can effect qualitative 
changes. From Section 3 onward we focus on the case of periodic orbits with a 
short region of sliding: this occurs close to so-called grazing- sliding bifurcations. 
First, in Section 3 we present a numerical observation for a minimal (2 -|- 1)- 
dimensional example based on the Hopf normal form. The example shows 
that, depending on the geometry of the problem, sliding may persist or not. In 
Section 4 we formulate a (still rather crude) general result about the persistence 
or destruction of sliding near grazing. Section 5 presents the e-cxpansions of 
the local return maps for the minimal Hopf example. Finally in Section 6 
we draw some initial conclusions and speculate along which lines we hope to 
generalize our results beyond the minimal example. 
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2 Simple example: destabilization of pseudo-equilibria 



Let us start with the simplest possible example of a Filippov slow-fast sys- 
tem. This example demonstrates the mechanism by which stable singular 
perturbations can destroy stable sliding in Filippov systems. Consider the 
system 

x^-sign[ex + (l-e)y] (7) 
ey^x-y. (8) 

In the quasi-static limit e = system (7) and (8) has the pseudo-equilibrium 
X — y — which is exponentially stable (even attracting in finite time). We 
recall that by a pseudo-equilibrium we mean an equilibrium of (7) and (8) 
that lies on the switching surface (a pseudo-equilibrium is typically not an 
equilibrium of any of the two flows near the surface). Consider the dynamics 
of the system (7) and (8) in the neighborhood of this pseudo-equilibrium for 
non-zero e. 



Let us express (7) and (8) in the general form given by (1), where 



/- 
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[yj 


_l{x-y)_ 



lx\_ 


-1 




\{.x-y)_ 



and the function h{x, y) that defines the switching manifold l-Lg is given by 

h{x, y) = 9x + {I - e)y. 

The vector normal to "Hs, is dh = [9^ 1—9]. The conditions for the existence 
of an attracting sliding region read (see [4]) 



dhf- > and dhf+ < 0, 
which gives for our example 

1-9 



Q 

9 H {x -y) >0, 



-9 + 



X — y) < 0. 



(9) 



(10) 



A subset of the switching manifold such that the vector fields /_ and /+ both 
point towards the switching surface along this subset (see Figure 1(b)) is 
called an attracting sliding region. A repelling sliding region is a subset of the 
switching manifold such that the vector fields /_ and /+ both point away from 
the switching surface along this subset. 

For 9^1 the switching manifold Hs can be parameterized by x, thus, the 
condition (10) simplifies to —e9 < x < e9 on Tig. Therefore, if 6' > then (9) 
holds along a segment of the switching line T-Lg in ^ neighborhood of the origin 
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Fig. 2. Two cases for system (7), (8). The pseudo-equilibrium either persists (case 
(a): 9 > 0), located on the sliding manifold Us = {Qx + (1 — 6)y = 0}, together 
with a small region where the system slides, or (case (b): < 0) it gets perturbed 
into a small periodic orbit and the sliding region near the origin is repelling. The 
manifold = {y = x} \s the slow manifold, on which the dynamics for e = 
lives. The manifolds M.± = {y = x it e} are perturbations of Mq. Each is invariant 
with respect to the ODE on its side of the switching manifold, containing no rapidly 
decaying dynamics. 

(see Figure 2(a)) and the stable pseudo-equilibrium persists for positive e and 
stays aX X = y = Moreover, in a sufficiently small neighborhood of the origin 
the trajectories converge to this pseudo-equilibrium following the sliding flow. 

If ^ < then (9) does not hold in the neighborhood of the origin and there is 
no stable pseudo-equilibrium (there is a repelling sliding region near the origin) 
and instead we have a symmetric exponentially stable periodic orbit around 
(0, 0) switching back and forth between /+ and /_. The point 

(■^switch; I/switch) (^(^1 ^ ^ Z/switch) 2/switch^ 

where the periodic orbit switches from /_ to /+ is parametrized by (see also 
Figure 2(b)): 

^_ 1 r + 2e-^ + e-^r-2 _ ^ e'^-l 

2 1 + e-^ ' y^-^'-^-^\-T^i 

where eT is the traveling time that the orbit spends on each side. Thus, for \9\ 
of order 1 (6^ < 0) the amplitude and the period of the periodic orbit are of 
order e; see Figure 2(b) for the qualitative picture. Consequently, the dynamics 
of the perturbed system (e > 0) are close to the dynamics of the unperturbed 
system (5 = 0) only in a very weak sense compared to the classical results of 
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[6]: instead of a stable fixed point we have an exponentially stable periodic 
orbit of small amplitude. This chattering phenomenon is similar to what is 
observed when the switching occurs with a small delay [2,7,12,13]. Example 
(7), (8) shows that chatter can also be introduced by fast dynamics even if the 
switching is perfectly instantaneous. 

The example also shows that singular perturbations may alter the defining 
quantities (9) of the sliding region by order 1 uniformly for s ^ and, thus, also 
change the dynamics qualitatively (in this example from exponentially stable 
pseudo-equilibrium to exponentially stable periodic orbit without sliding) . This 
will be further highlighted in Section 5 where a periodic orbit in a singularly 
perturbed Filippov system is analyzed in detail. 

Figure 2 shows the arrangement of invariant manifolds and the switching 
manifold for the example (7), (8). For each of the two flows there exists an 
invariant manifold M.± consisting of all trajectories (in the 2D-example just 
one trajectory) without rapidly decaying part. These manifolds Ai± are both 
0(e)-pcrturbations of the so-called slow manifold Aio, defined by setting e = 
in the equation (8) for the fast dynamics. In the example the manifolds are 
A4± = {y = X ± s} and A4q — {x = y}. Every trajectory that spends a time 
of order 1 in the subdomain — {{x, y) : h{x, y) > 0} ends up exponentially 
close to the invariant manifold A^+: its end point has a distance of order 
exp(— c/e) from (similarly for M._ in = {{x,y) : h{x,y) < 0}). The 
example shows that the invariant manifolds Ai+ and Ai- typically differ by 
a term of order e. This implies that any trajectory that switches back and 
forth between T-L^ and T-L- shows a small boundary layer immediately after 
switching: after switching the trajectory has a short time interval during which 
it relaxes to the invariant manifold of the other flow. These small boundary 
layers occur also after switching to sliding and have been studied already in 
[8,9,10]. In example (7) for ^ < both parts of the periodic orbit are still in 
the boundary layer. 



3 Grazing-sliding in singularly perturbed Filippov systems — il- 
lustration 



In order to understand the effect of stable singular perturbations on periodic 
orbits with sliding we focus on periodic orbits with a short sliding segment. 

One common scenario for periodic orbits with short sliding segments is the so- 
called grazing-sliding event; see [4] for a classification of discontinuity-induced 
bifurcations. This phenomenon is a codimension-one event, which can be 
observed generically if the system depends on one additional parameter. 
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n. 




(a): /i < 1 



(b): ^i=l 




(c): II > 1 



Fig. 3. Grazing-sliding scenario for minimal example (11). (a): the periodic orbit lies 
entirely in 71+; (b): the periodic orbit touches the switching line {xi = —1} at Xg] 
(c): the periodic orbit has a small sliding segment from Xh to xi. The line Ci is the 
Poincare section chosen for the return map Mi in (12). 

3. 1 Grazing-sliding — illustration for minimal example 



A minimal example for the grazing-sliding event, the Hopf normal form for /+ 
combined with a constant vector field for /_, is 



X = < 



where we keep u ^ fixed and vary fi (which is always greater than 0) from 
below 1 to above 1 (see Figure 3). The switching line for this example is 

h°(x) = Xi + 1, 

and the flow in 7/+ = {h^ > 0} has a unique stable periodic orbit as shown in 
Figure 3(a). In fact, the flow in 7/+ can be decomposed into a pair of uncoupled 
equations for the polar coordinates r and of a; G {x = r ■ (cos*/?, simp)): 

r = fir — r^, (p = 

and the periodic orbit is a circle around the origin corresponding to r = y^, 
if = ujt. This periodic orbit is also a periodic orbit of the Filippov system (11) 
for yU < 1. When we change /i to values larger than 1 the periodic orbit of 
the Filippov system (11) acquires a sliding segment starting from some point 
Xh = [—l,Xh,2Y {^2,h > 0) and ending at some point xi = [— l,x/,2]"^ (a^«,2 < 
is defined by the condition dhf_^{xi) = 0), making the overall orbit piecewise 
smooth (shown in Figure 3(c)): it has a corner at Xh and its second derivative 
is discontinuous at xi. The periodic orbit is grazing the switching manifold 



fXXi — LOX2 
(jJXi + jJ,X2 
1 





if xi > -1 



if xi < —1, 



(11) 
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H'^ — {h^ = 0} ior fjL — 1 as shown in Figure 3(b) in a point Xg, which satisfies 
several conditions simultaneously: 



(1) The point Xg lies on the stable periodic orbit of the flow in and on 
the switching manifold H^. 

(2) The flow in T-t+ is quadratically tangent to the switching manifold {h — 
and h > in Xg along trajectories following /+), and 

(3) the flow in points toward the switching manifold in Xg. 

The dynamics near the periodic orbit for ?a 1 is described completely by its 

local return map (also called Poincare map). This map is defined by following 
the fiow from a small cross-section transversal to the periodic orbit back to itself. 
We choose the local return map to the section Ci = {x : xi = 0, X2 > 0} 
(see Figure 3(b)), which has for // = 1 the form 

(^-4nM _ 1] + 0(\X2 - IP) if X2 < 1 

M,X2 = l+\ ^ ^' ' ' ^ \ (12) 

10 if X2 > 1. 

Thus, the map Mi is picccwisc asymptotically linear in its fixed point X2 = 1 
corresponding to the periodic orbit. For /i > 1 the map is constant near 
its fixed point. In summary, in system (11) the only effect of the grazing under 
variation of /i is that the periodic orbit changes its shape, and that its Floquet 
multiplier jumps from exp{—47rfi/uj) to at the grazing parameter /j, — 1. The 
periodic orbit is stable for all ~ 1. 



3.2 Singular perturbation of grazing-sliding — observations for minimal ex- 
ample 



Now let us consider a stable singular perturbation of system (11). We will 
couple (11) with a one-dimensional stable fast subsystem for a fast variable y 
where the coupling also occurs in the switching function. Thus we obtain 



X 



ey^e 



fiXi — UIX2 — {x\ + x\)xi 

UOXi + 11X2 — {x\ + X\)X2 

1 



jlXi — UJX2 — {x\ + x\)xi 



if Oxi + (1 - e)y > -1 

if exi + {l-e)y< -1 
+ xi-y. 



(13) 
(14) 



The quasi-static approximation replaces y by yQ{x,fi) = Xi (e = in (14)). 
Thus, the reduced model is (11) which is identical with the slow part (13) of 
the singularly perturbed system except that the switching function depends 
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(a) e = -0.5 (b) e = 0.5 



Fig. 4. Approximate local return map to £i, P, for e = 0.01, ^ = 1, w = Svr obtained 
numerically as the first return to Ci (neglecting the distance of the return point to 
the slow invariant manifold = {(x, y); x\ = y} of the flow in %+)■ Intersections 
with the diagonal indicate fixed points corresponding to period-one orbits. 

now in part on the fast variable y ioi 9 ^ 1. Namely 

h{x, y, e) = 6xi + (1 - 6)y + 1. 

The question is: how does the perturbation from the reduced system (11) to 
the slow-fast system (13), (14) affect the local return map Mi of the periodic 
orbit at /i = 1? For /i = 1 the grazing periodic orbit drawn in Figure 3(b) exists 
also in the slow-fast system: its grazing point is now pg = (—1,0,-1), and 
on the orbit the fast variable is given by y{t) = xi{t). We make the following 
observation (see also Figure 4): 

Observation 1 (Perturbation of local return map near grazing) 

The perturbed return map for e > is a one- dimensional return map P to the 
line Ci = {{x, y); Xi = y = 1} if we ignore terms of order exp(— c/e). 

If 9 < then the map P is discontinuous for e > 0, and the discontinuity is of 
size 0{e). 

If 6 > then the return map P is continuous also for e > but the slope near 
the point Pg is of order 1 uniformly for e — >■ and X2 > 1- 

Figure 4 shows how the return map Mi of the grazing periodic orbit is perturbed 
for positive e, and 9 <0 (Figure 4(a)) and ^ > (Figure 4(b)). The grazing 
point Pg corresponds to X2 = = 1. For X2 < 1 the return map P is identical 
to Ml. We observe that for ^ < the return map becomes discontinuous at 
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X2 — 1 and that the size of the jump is of order e. Also, the hmiting slope 
for X2 > 1 appears to be — oo. For > wc note that the limiting slope for 
^ +1 is not of order e. The slope is negative and of order 1. For comparison, 
the map Mi is identical to 1 for a;2 > 1. The change of shape of the graph 
of the return map from Mi to P has consequences for the dynamics near the 
grazing orbit. For ^ < the period-one orbit disappears as fj, passes through 1 
from below and we observe an inverted period-adding cascade [5] . For ^ > 
the period-one orbit persists but may change its stability depending on the 
slope of the graph in Figure 4(b) for X2 — )■ +1. The local attractor (if any) 
for // > 1 also depends on this limiting slope (see [4] for a recent textbook 
overview of possible scenarios in piecewise locally linear maps). 



4 SingulEir perturbation of greizing-sliding events 

This section gives a criterion that helps to distinguish whether a given periodic 
orbit that grazes the switching surface in a general slow-fast Filippov system 
undergoes a bifurcation as shown in Figure 4(a) (that is, the local return map 
is discontinuous) or if we have a grazing-sliding case shown in Figure 4(b) 
(the local return map is continuous). Let us assume that the reduced Filippov 
system (6) depends on a parameter /i, and that (6) has a family of periodic 
orbits x{t]ix). We assume without loss of generality that the periodic orbit 
x(t; n) has no sliding segment near t = for /x < 0, and it lies entirely in the 
region "H^ = {x : h(x,yo{x),0) > 0} locally near t = (see Figure 3(a)). The 
precise characterization for a grazing-sliding event for the family x(t; fi) at a 
parameter value /xq (/xq — without loss of generality) can be given by using 
the function xo '■ (^;A*) ^ h{x{t; iJ,),yo(x{t; iJ,)),0). The function xo gives the 
value of h depending on time along the periodic orbit. The family x{t; ji) grazes 
at /X = if there exists a time to near such that 

(1) Xo(^0;0) = (x(io;0) is on the switching manifold Hs)-, 

(2) 9iXo(^0;0) = (the orbit grazes l-Ls in ^o), 

(3) 9iXo(^OjO) > (the grazing is quadratic), 

(4) (92Xo(^0)0) < (the parameter /i unfolds the grazing), 

(5) [dih + d2hd2g^^dig\f- > in (a;(to; 0), yo(a;(to; 0)), 0) (the reduced vector 
field in %_ = {x : h{x,yQ{x),0) < 0} points toward switching manifold). 

(The functions h, f± and yo may all depend on /i in this list.) Note that this is 
a codimension-one event because we have two equality conditions (1) and (2) 
and can adjust the time tg- The conditions 1-5 guarantee that the periodic 
orbit and the switching manifold interact as shown in Figure 3(b). See also [4] 
for a detailed discussion of grazing-sliding events for periodic orbits. 
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Purthcrmorc, let us assume that the periodic orbit x{t; n) is entirely in Ti^ 
for fi < (not only locally near to)- Then the standard singular perturbation 
theory of [6] guarantees that for sufficiently small e > the full system (3), (4) 
also has a family of periodic orbits {x, y){t] ii, e) which depends smoothly on /i 
up to some parameter value — 0{e) where it grazes quadratically at a point 
(xs, He) = {x, y){to', fJ'e, We Can keep the time point of the grazing equal to to 
independent of e without loss of generality by shifting the time in the family of 
periodic orbits because the slow-fast system (3), (4) is autonomous. This means 
that the function Xe ■ {t-ilj) ^ h{x{t; ii,e),y{t; iJ,,e),e) also satisfies the first 
four conditions 1-4 for /J, — and ^q. This is a simple persistence argument 
because only the smooth vector field /+ is involved. However, in general the 
equivalent of condition 5 (that the vector field on the other side points toward 
the switching manifold) is not necessarily satisfied for the full system (3), (4) 
even though it is satisfied for the reduced system. 

The following lemma states how this affects the local return map of the periodic 
orbit at /i^: 

Lemma 2 (Grazing-sliding of periodic orbits) Suppose that the reduced 

system (6) depends on a parameter fi and has a grazing-sliding periodic orbit, 
lying entirely in and satisfying the conditions 1-5 for the parameter value 
11 = at xq = x{tQ; 0). 

Then for sufficiently small e > the family of periodic orbits in the full system 
(3), (4) also grazes for a nearby parameter value /le — 0{e) and we have two 
generic cases. 

Case 1 If dihf- — dihf+ > in (xq, yo(xo), 0) then for sufficiently small 
e > the local return map of the grazing periodic orbit is continuous (and 
piecewise smooth). 

Case 2 If dihf- — dihf+ < in {xo,yo{xo),0) then for sufficiently small 
e > the local return map of the grazing periodic orbit is discontinuous in 
the fixed point corresponding to the grazing periodic orbit. 

In our set-up the local return map does not depend on the choice of the cross- 
section as long as the cross-section is not taken at to (return maps to different 
cross-sections are equivalent and the coordinate change is a diffeomorphism 
obtained by following the flow in 7^+). 

As the minimal example (6) shows, a typical feature of the grazing-sliding event 
is that the local return map along the periodic orbit is piecewise smooth in its 
fixed point corresponding to the periodic orbit x{-; /i) for the parameter n = lie 
at which grazing occurs. More precisely, the return map is smooth (with a 
uniformly bounded derivative) everywhere except along a manifold containing 
the fixed point. At this manifold the map is only (Lipschitz) continuous. The 
lemma states that in the second case the return map of the grazing orbit in 
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the singularly perturbed system develops a discontinuity for positive e at this 
manifold. The appearance of a discontinuity is guaranteed because we assume 
that the grazing orbit lies entirely in one domain, always following the flow 
(/+, g) (see Appendix A for details). In Section 5 we show that for the minimal 
example (13), (14) the size of the discontinuity can be expected to be of order 
e. 



The quantity distinguishing case 1 and case 2 can be evaluated without knowl- 
edge of the right-hand-side g of the fast dynamical subsystem (4). This means 
that one can check if sliding persists under stable singular perturbations of a 
term Uq^x) in the reduced model without having to know how this perturbation 
looks like. 



The two cases are distinguished by checking if the vector field in the domain 
%- = {{x,y) : h{x,y,£) < 0} points toward the switching surface T-Lg near 
the grazing point. It turns out that in the second case the reason for the 
discontinuity of the return map is the presence of a small repelling sliding 
region. For the first case the sliding condition is satisfied such that sliding 
occurs and the general theory developed in [4] (Theorem 8.1) applies. The proof 
of Lemma 2, given in Appendix A, simply has to check this sliding condition. 
Theorem 8.1 of [4] can also be used to approximate the one-sided derivatives 
of the local return map of the full system (3), (4) for the sliding case (case 1) 
of Lemma 2. We show in Section 5 for the minimal example (13), (14) that 
the limit for £ — > for these one-sided derivatives can be different from the 
one-sided derivatives obtained for the reduced model. 



Lemma 2 is rather crude: a more detailed analysis is required to find out how 
the return map of the grazing periodic orbit actually looks like and how its 
features (for example, the size of the discontinuity or the one-sided derivatives) 
depend on e. This analysis is technical for a general system such as (3), (4), 
in particular if the dimension of the fast variable y is larger than 1. Thus, we 
construct the ^-expansions for the return map only for the minimal example 
(13), (14). The example allows us to address both cases of Lemma 2 by varying 
its system parameter 6. As long as the fast variable is one-dimensional, the 
generalization to arbitrary singularly perturbed systems with a grazing-sliding 
periodic orbit is straightforward. 
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Fig. 5. Illustration of the manifolds near the grazing point pg for the slow-fast 
example (15), (16) for e > 0, 6 < and = 1 (zoom-in near the grazing point pg 
of the periodic orbit). The grazing periodic orbit lies entirely in T-L+ (and, thus, in 
The other invariant manifold Ai- (for T-L-) has a distance of 0(e) from A^+. 
The illustration shows the grazing periodic orbit and a typical (for 9 < 0) switching 
trajectory near the periodic orbit, switching from to back to E^, and then 
re- approaching A^-|_. 

5 Expansion of return maps for the minimal example 



Let us return to the minimal example 



X 



jlXi — (jJX2 — {x\ + x\)Xi 
OOXi + /iX2 - {x\ + X2)X2 



if exi + {l-d)y> -1 
if Oxi + {l-e)y < -1 



llXi — UX2 — {xl + x\)xi 



+ xi-y 



(15) 
(16) 



where a; > 0, < 1 and /i ~ 1. We have chosen the right-hand-side g of the 
fast equation (16) such that for e > the subspace = {{x^y) : xi = y} is 
the exact slow invariant manifold of the flow E"^ (following x = fj^, ey = g). 
That is, any trajectory will converge to with a convergence rate of order 
as long as it stays in the half space 7/+ = {{x, y) : Oxi + (1 — 0)y > — 1}. 
The graph of the invariant manifold Ai+ for all and e is given by 



Xl. 



(17) 



Since /+ does not depend on y the stable fibres of the manifold Ai+ are: 
J^+{xq) = {{x,y) : X = Xq, y ^ M} for all Xq G (all points within a stable 
fibre converge to each other with a rate 0{e~^) following fiow forward in 
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time). The graph of the slow invariant manifold M.- oi (the flow following 
X — f^, ey — g) is not known analytically. Its expansion up to order e is 



yjn-{x, = Xi + £ fXXi - UX2 - {Xi+ X2)Xi - 1 +0{6). (18) 

The expressions for the slow invariant manifolds for and E-, (17) and 
(18) differ from each other by a term of order e; the difference is e + 0{e'^) 
for X = (—1,0), y = —1, 11 = 1, thus, it is non-zero at the grazing point 
Pg = (—1,0,-1) at /X = 1. Hence, we expect that any trajectory crossing 
the switching manifold will show a small boundary layer; see Figure 5 for 
an illustration of the invariant manifolds A4± and the switching manifold 
Us — {{-c, y) : 9xi + (1 — 0)y = —1} near the grazing point at = 1. 

The switching manifold Hs for the example (15), (16) is a two-dimensional 
plane. Its intersection with the invariant manifold is 

nsnM+ = {{x,y):xi = y = -1}. 



The two cases of Lemma 2 correspond to ^ > (Case 1) and ^ < (Case 2). 
They have to be treated differently and result in qualitatively different return 
maps as Figure 4 shows. The return maps for both cases are superpositions of 
a smooth global return map (following E+) and a correction in the vicinity 
of the grazing, the Poincare-section discontinuity mapping (PDM, [4]). In the 
following we derive an expression for the return map to the section Si = 
{{x,y) : = 0, a;2 > 0, 7/ e M} which is away from the grazing point (see 
Figure 5 for £1 = Ei fl M.+). More specifically we consider the local return 
map to a smaU neighborhood W C Si of the point (x, y)^ = (0, 1, 0)^, the fixed 
point at /i = 1 corresponding to the grazing periodic orbit. Any trajectory 
starting in U spends a time of approximately 7r/(2w) following E^ before it 
reaches the vicinity of the switching manifold Hs in a small neighborhood V 
of the point Pg = (— 1,0,— 1)^. During this time the difference between y and 
Xi (which is the distance to the slow invariant manifold A4+) decays such that 
|y — xi| ~ exp(— 7r/(2c(;e)) when the trajectory reaches V, which is beyond all 
orders of 5. A trajectory leaving V follows E^^^ for a time of approximately 
371 /{2u) until it reaches Si. Again, after this time the y-component and the 
Xi-component of the trajectory will be exp(— 27r/(3a;£:)) close to each other. 
This means that the return map P is a one-dimensional map from the line 
>Ci = Si n Ai+ = {{x , y) : Xi = y = , X2 > 0} back to itself if we ignore terms 
of order cxp{—c/e) where c > is of order 1 (the return line Ci is also indicated 
by a dashed line in Figure 5). The overall return map P is a composition of 
four maps. CaUing the time derivative of h with respect to each of the flows as 
h'_^, respectively, 

h'^{x,y) = dh{x,y)^El{x,y)\t=o = [dihf± + e~^d2hg]{x,y), 
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(dropping the argument e) the maps Pi to P4 act as follows 

(1) Pi, maps from the line £1 = {xi = y = 0, X2} HU to the curve £2 = {xi ~ 
y, h'_^{x,y) = 0} n V by following the flow within the slow invariant 
manifold M.+ (see Fig. 5), 

(2) P2, the Poincare-section discontinuity mapping from £2 to the surface 
S = {{x,y) : h'_^{x,y) = 0} fl V, which contains the curve £2, 

(3) P3, maps from E back to £2, following the projection along the stable 
fibres of P3{xi,X2,y) = {xi,X2,xi), 

(4) P4, maps from £2 C V back to £1 C W following inside the slow 
invariant manifold 

The composition 

P = P4 o P3 o P2 o Pi 

maps the line £1 back onto itself and it is a one-dimensional approximation 
of the true two-dimensional return map near the grazing periodic orbit up to 
terms of order exp(— c/s). The map P is conjugate (up to the diffeomorphism 
Pi) to the map 

[Pi o P4] o [P3 o P2] =: Pgiob o Pdm, 

which maps from £2 back to itself. The smooth map Pgiob = Pi o P4 : £2 1— )■ £2 
is the global return map around the periodic orbit following the smooth flow 
E^. Thus, Pgiob has the fixed point 

2;i,giob = -Va* (19) 

corresponding to the stable periodic orbit of E+. Its hnearization at the fixed 
point is stable and is identical to the non-trivial Floquet multiplier of the 
periodic orbit of which is exp(— 47ryu/w). This means that the smooth 
global return map Pgiob = Pi o P4 is described to first order by 

-Pgiob(a;i) = -v^ + exp {-Ann/u) [xi - xi^giob] + 0{[xi - xi^giob]^)- (20) 

The other part of the return map P, the projected discontinuity mapping 

Pdm = P3 o P2 

depends on the sign of the parameter 9 in the switching function h. For each sign 
of 6 we study a range of fi, x and y near the grazing point (x, y) = Pg = (—1, 0), 
/i = 1. The size of the range depends on e: it is 0{e) for 6* < 0, and it is 
0(£^) for ^ > 0. We take this into account by coordinates corresponding to 
a zoom-in into the neighborhood of the the grazing point Pg (blowing up the 
small neighborhood to size 1). 

Lemma 3 (Discontinuity mappings) 
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Case 6 < 0: Introducing the coordinates and q by xi — — 1 + and 
II — 1 + eq, the map Pqm has the form 



-Pdm(Ci) 




+ 0{e) i/6<0. 

(21) 

Case 9 > 0: Introducing the coordinates and q by Xi — —1 + £^^i and 
II — 1 + e^q, the map Pdm has the form 



-Pdm(^iJ 



6 «/6>o 



(22) 



where the graph s{6,uj, ^i,e)^i is parametrically defined by the sliding time 
is > 0. The zero-order term can be represented as 



s{e,uj,iumi 



CO 



2 

":i- 



(1 - ^) [1 - eiq>{-ts/0)] + ts 



9 + {l- 

e-[e + t, 



9)exp{-ts/e) 

\exp{-tJ9) 



(23) 



^ + (1 - ^) ex.p{-ts/e) 



The quantity So{9) appearing in (21) is a uniformly positive constant (only 
depending on 9). It is implicitly defined by the equation (B.15) in Appendix B.l 
(see also graph of So(^) in Figure B.2). The remainder term 0{e) in (21) for 
^1 < may contain small corrections to the constant and square-root terms in 
^1. However, the form (21) guarantees that for sufficiently small values of the 
singular perturbation parameter e the overall return map Pgiob ° -Pdm changes 
from a piecewise linear map for the reduced model (e = 0) to a discontinuous 
map for the full model {e > 0). The size of the jump at the discontinuity is 
of order e in the original coordinates. Moreover, the slope of the map next to 
the discontinuity is infinity from one side. Figure 4(b) shows a numerically 
computed graph of the approximate overall return map P to jCi for ^ = —0.5, 
e = 10~^ and ii = I near the critical value of the coordinate X2. This graph 
neglects exponentially small terms (the numerical distance of the return value 
to the slow invariant manifold is of the order 10~^° at the return section). 
Appendix B.l gives the details of the derivation of the expansion of the map 
Pdm in (21). 

In (22) the function s is smooth in all arguments such that higher order terms 
do not destroy the continuity of the map Pdm- Figure 6 shows the graph 
of s(0.5,87r,^i,0). The value of s at = (that is, t^ = 0) is {9 - l)/9. 
Appendix B.2 gives the details of the derivation for the discontinuity mapping 
Pdm in (22). The change of the slope s from its value at = to slope for 
^1 — )■ — oo (as predicted by the reduced model) occurs over a range of order e"^ 
in the original coordinates after grazing. Figure 4(a) shows the approximate 
overall return map P to Ci for 9 — 0.5, e — 10~^ and n — 1 near the critical 
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-0.99 




-1 



-0.03 -0.02 -0.01 

Fig. 6. Illustration of the implicitly defined graph s{9,uj,-,0), as defined by (23) 
(parametrized by tg, to = Sw, 9 = 0.5). 

value of the coordinate X2- 



6 Conclusion 

Stable singular perturbations have a much stronger influence on the dynamics 
in Filippov systems than in smooth dynamical systems. We demonstrate 
that stable pseudo-equilibria and stable periodic orbits with sliding do not 
necessarily persist. We study periodic orbits with an infinitesimally small 
sliding segment, that is, close to a grazing-sliding bifurcation. We found two 
generic cases depending on the geometry: the local return map around the 
grazing periodic orbit develops a discontinuity if the condition on the existence 
of an attracting sliding region is violated. Otherwise, the continuity of the 
return map persists but the asymptotic slope may have a change of order 1 
(uniform for s ^ 0). 

The qualitative change of the local return map induces qualitative changes to 
the dynamics on a small scale. A piecewisc discontinuous map with a square- 
root singularity of the slope on one side of the discontinuity, as occurs for ^ < 
in the minimal example, shows inverted period-adding cascades of periodic 
orbits if one varies the parameter /i through its critical value [5] . The parameter 
range where these cascades can be observed is of order e. In the other case, 
^ > 0, the observed dynamics depends strongly on the one-sided derivative 
s{6, 0,0) defined in Lemma 3. It can be chaotic if s{9, 0,0) < —1, which is 
possible for small 9. Our analysis is valid on a scale of order e"^ in phase space 
and parameter space. 

The results of Lemma 3 can be generalized to higher-dimensional slow-fast 
systems in a straightforward manner as long as the dimension of the fast 
subsystem is 1. There are, however, some technical difficulties to generalizing 
the expressions of Lemma 3 for higher dimensional fast subsystems (y e R"*, 
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m > 1); a trajectory following the dynamics inside the stable fibres (following 
a linear stable ODE) may intersect the switching hyperplane several times. For 
the expressions (21) and (22) in Lemma 3 we exploited that this is not the 
case. In R-*- every trajectory in a stable linear system approaches the origin in 
a monotone (increasing or decreasing) fashion, which is not true in in the 
Euclidean norm. Furthermore, the Poincare-section discontinuity mapping is 
only implicitly given as the root of a nonlinear equation even for the minimal 
example. In general, this implicit expression is determined by the intersection 
of a trajectory following a stable linear system with a hyperplane. 

What are the consequences of the small-scale instabilities described in our 
paper in practical applications? They offer one explanation for discrepancies 
between predictions of a simple model and real observations: if the model is a 
smooth dynamical system and predicts stable periodic motion then observed 
noise around the predicted periodic motion in real observations is evidence for 
the presence of (unmodelled) inhomogeneous forcing. Fenichel's Theorem [6] 
guarantees that unmodelled strongly stable degrees of freedom (that is, the 
fast subsystem) cannot be the cause of this apparent noise. If the model is a 
Filippov system and sliding occurs in the predicted periodic motion then the 
fast subsystem can disturb the periodic motion in a way that it is no longer 
periodic (for example, chaotic, or periodic with a much higher period). Hence, 
if the model predicts sliding motion then apparent noise in the observation is 
not evidence for unmodelled inhomogeneous forcing. 

The influence of fast dynamics also provides a mechanism for chatter (that 
is, rapid sequences of switches instead of sliding) which competes with other 
possible mechanisms. For example, switches in electronic circuits are known 
to be non-ideal [1] . Suppose that an electronic switch acts with a small delay 
T but is otherwise ideal. Taking this delay into account one would predict 
similar effects as described in this paper [7,12,13]. Similarly, if the switches 
are implemented by fast actuators (with speed 1/r as described in [9]) one 
would expect to observe chattering with frequency of the order 1/r. However, 
if unmodelled capacities (or other stable degrees of freedom) give rise to fast 
stable subsystems with decay times e > t then the frequency of the chatter 
will be limited by 1/e, which is smaller than 1/r. 

Another open question is if mechanical systems with dry friction have a special 
structure that prevents singular perturbations from destroying sliding (the case 
^ < in our examples). Singular perturbations in mechanical systems may 
certainly change the stability of sliding periodic orbits near grazing. This is also 
possible in the case corresponding to ^ > in our examples: sliding persists 
but at grazing the one-sided derivative of the grazing orbit has a change of 
order 1 due to the fast dynamics. We demonstrate this effect in Appendix C 
for the classical autonomous dry-friction oscillator. 
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A Proof of Lemma 2 



The assumptions of Lemma 2 imply that the system 



has a family of periodic orbits (x(t; iJ,,£),y{t; ji^e)) that grazes the manifold 
T^ti = {h = ^} for all sufficiently small e > 0: this follows from the fact that the 
reduced vector field [e = 0) has a family of periodic orbits {x{t; fi, 0),y{t; fi, 0)) 
which grazes at /i = 0, and that the slow- fast system (A.l) and (A. 2) on 
its slow invariant manifold in the domain T-L^ = {/i > 0} is a regular 
perturbation of the reduced vector field. Without loss of generality we can 
(by reparametrizing /i) assume that the grazing occurs ai /i = for all e, that 
the time at which grazing occurs on the orbit is t = for all e > 0, and 
that the point (x{0,0,e),y{0,0,e)) in which the orbit grazes is {x,y) = (0,0) 
(by shifting time of the autonomous system and shifting the origin of the 
coordinate system). In our new coordinates the family {x{t] /i, s),y{t] n, e)) is a 
family of periodic orbits, which depends smoothly on both parameters (//, e) in 
the parameter interval [— /UmaxiO] x [O.emax] (both, /Umaxj and Smax arc chosen 
sufficiently small). The periodic orbits are entirely inside and graze T-Ls 
quadratically in (0, 0) for = (for all s e [0,£max])- 

This quadratic grazing condition means that the function 



(which is smooth on R x [— /Xmax, 0] x [0, £max]) satisfies for all e e [0, £max]: 

(1) x(0, 0, e) = (Tis contains the origin for n = 0) 

(2) 9ix(0, 0, e) = (the periodic orbit grazes Tig in the origin for /jl — 0) 

(3) 9^x(0, 0, £) > (the grazing is quadratic) 

(4) 92x(0, 0, £) < (/X unfolds the grazing). 

In order to find a local return map of the periodic orbit {x{t,0,s),y{t,0,s)) 
we have to choose an appropriate Poincare cross section: the conditions 2 and 




\f_{x,y,fj.,e) if h{x,y,fj,,e) <0, 
\f+{x, y, fi, e) if h{x, y, n,e)>0 
g{x,y,ii,e) 



(A.1) 



ey 



(A.2) 



X ■ it,ii,e) ^ h{x{t;ii,e),y{t;ii,s),ii,s) 
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3 guarantee that the set 



S := {{x,y) : ed^hU + d2hg = 0, \\{x,y)\\ < 6} 

defines a smooth hypersurface near the origin, which intersects the periodic 
orbit and the switching manifold Tis transversally for all sufficiently small e 
and /i. The set S is the set of local minima of h along trajectories of the flow 
£"+ generated by the vector field {f^,g/s). We define our Poincare section as 
El := where T'(//, s) is the period of the periodic orbit and is 

the time-i map of the flow E^. 

Theorem 8.1 in textbook [4] states that, if the vector field [f^^g/e) points 
toward ?^s(0,e) in the origin then the local return map to Ei is continuous 
and it can be written in the form 

P = E^-Si ^ p^^ ^ ^T(M,.)/2 

where E^^^^ is the map obtained by following from E to its first intersection 
with El (this is locally a diffeomorphism, the traveling time is close to T(/i, £)/2). 
The map Pdm maps E back to itself and accounts for all effects due to the 
discontinuity (it is called the Poincare discontinuity map in [4]). Whether the 
vector field {f-,g/e) points toward or away from the switching manifold in the 
origin is determined by the sign of the function 

■ {x,y,l^,e) ^ [dihf-+e-^d2hg]{x,y,i^,e) (A.3) 

in the point {x,y,fi,e) = (0,0,0,6). Due to condition (2) we know that the 
time derivative of h when following E^ vanishes in (0, 0, 0, e): 

dihU + £-^d2hg ^0. (A.4) 

Thus, combining (A.3) and (A.4) we obtain that the quantity distinguishing 
the two cases in Lemma 2 is in fact 



(7_(0,0,0,£) = - dihU. (A.5) 

That is, in case 1 the vector field E^ points toward the switching manifold Hs 
and we can apply Theorem 8.1 from [4], proving the claim of Lemma 2 for 
case 1, whereas, in case 2, E_ points away from the switching manifold Us- 
Theorem 8.1 from [4] also gives an asymptotic expression for the one-sided 
derivatives of the map Pdm in (0,0). The example from Section 5 will show 
that the limit of these one-sided derivatives for £ — >■ can be different from 
the one-sided derivative when one apphes Theorem 8.1 from [4] to the reduced 
system {e — 0). 

It remains to be shown that in case 2, when the vector field {f_,g/e) points 
away from T-Ls in the origin for fJ, — 0, the return map to Ei is discontinuous 
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for £ > 0. Let us define the sets 

H+ = Ef^^''^^'{{x,y) e S : hix,y,e) > s} 
N = El^^^{{x,y)ei::h{x,y,6)<0}. 

The set Hi consists of all points in the Poincare section Si for which Pbm is 
the identity (that is, trajectories starting in Hi do not switch to the flow 
during the evaluation of the return map P). The set consists of points in Si 
that do not have a pre-image (they are images of a repelling sliding section 
under Ef^^'). 

The set PHi C Si can be written as 

PH+ = y) e S : h{x, y, 0) > 0} 

since E^'^^'^ o ^^('^''^^/^ jg the whole return map P when Pdm is the identity. 
Thus, the set PHiUN is an image of the whole cross-section S under the local 
diffeomorphism . Hence, PHi U N contains an open neighborhood of 

the fixed point po e Si corresponding to the grazing periodic orbit. Moreover, 
P is a diffeomorphism between Hi and PHf . This implies that points from 
Si \ Hf cannot be mapped into the neighborhood of po under P (they can 
neither be mapped into N nor into PHf). However, po is on the boundary of 
Hi and, thus, it is in the closure of Si \ Hi . Consequently, there exist points 
arbitrarily close to po that do not get mapped into the neighborhood of po 
under P, which proves that P is discontinuous in case 2. 

We note that the proof of discontinuity in case 2 of Lemma 2 is not specific to 
slow-fast systems (we only use that the vector field E_ points away from the 
switching manifold at the grazing point (0, 0, 0, e)) but it relies on the fact that 
the grazing orbit lies entirely in the domain governed by the flow P+ (that is 
why P restricted to Hf is a diffeomorphism). In general, for example, if the 
grazing periodic orbit has sliding segments of positive length in other parts of 
the phase space, the return map may not be discontinuous in case 2. 



B The Poinccire discontinuity mappings of the minimal example in 
Section 5 

This appendix derives the dominant terms of the projected discontinuity 
mapping 

-Pdm = Pi o P2 

as defined for the minimal example in Section 5. The map Pdm maps the line 
'^2 = {h'j^{x, y) = 0, Xi = y} back to itself in the neighborhood of the grazing 
point Pg — (— 1,0,— 1). Orbits of crossing £2 in a point with coordinate 
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y 




Fig. B.l. The Poincare discontinuity map for the case < near the point pg (0 in 
the rescaled coordindate system (B.2)). 



Xi > — 1 do not switch to or slide: by definition of £2 eacli orbit has its 
local minimum of h exactly when it crosses £2, and h is identical to xi + 1 
on the invariant manifold M.+ of Ej^. Thus, £2 can be parametrized by the 
coordinate Xi, and for Xi > — 1 the map Pdm is the identity: Pt)m{xi) = Xi. 
For xi < — 1 the map Pdm is itself again a composition of several maps: 



• Hi follows the flow backward in time from £2 to the switching manifold 
Hs, specifically, the line £3 = {xi = ^ = 0}, 

• the map 112 depends on the sign of the parameter 6: 

6 < 0: II2 follows E_ from £3 forward in time until one hits the switching 

manifold Hs again; 
6 > 0: II2 follows the sliding flow Eg, defined by (2), starting from £3 and 

staying in the switching surface Tig until the flow E^ becomes tangent to 

Us, that is, h'_^{x,y) becomes zero; 

• P3 projects along the stable fibre of P+ down onto the slow invariant manifold 

M + : P3(Xi,X2,l/) = (Xi,X2,Xi). 

• 113 follows -E+ on (forward or backward in time) to £2. 

We have swapped P3 and 113 in (B.l), which is possible because P3 and II3 
commute. 



Pdm = lis o P3 o 112 o III 



(B.l) 



where 
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B.l The case 6 <0 



For ^ < the parameter range of interest is of order e. Similarly, the scale 
of the dynamics transversal to the limit cycle (that is, the xi and the y 
components) is of order e. We zoom into the neighborhood of the grazing 
point pg = ( — 1,0, —1) by introducing the new quantities ^, iJj and q (which 
are rescalings of the coordinates x, y and /x): 



xi^ -1 + 5^^1 



(B.2) 



Any point in the phase-parameter space with ^, ip and q of order 1 is in the 

vicinity of the grazing bifurcation. The scaling (B.2) squeezes the direction along 
the grazing periodic orbit by scaling it only hy 6 = y/e. In the new coordinates 
the switching manifold 'H,, is {(^,V') : + (1 - 9)4' = 0}. System (13), (14) 
in the new coordinates has the form (with rescaled time tnew = and 
truncating terms of order 5^): 



-a;6 + [ei-2ei-g]5 + 0((53) 

1/6 




if 9^1 + {l-9)ij>0 



if 9^1 + {l-9)ij <0 



5^ = -V' + - co^2S + - 2^1 - q]S' + 0(5^). 



(B.3) 



(B.4) 



We expand the quantities h'_^ and h'_ in the new coordinates in 6 (dropping 
0{S^) and substituting V by 9^i/{9 - 1)): 



5-'h'_ 



9 + ^1 + {9- l)uj^2S 
+ {l-9) [e,-q-2^,] 6' + 0{5'). 



(B.5) 



Thus, in the vicinity of the grazing point the vector field points toward the 
switching manifold Hg if /i+ < 0, that is, 

Ci<ujC2S-[e2-q]S' + 0{5'), (B.6) 

and E_ points toward T-Lg if h'_ > 0, that is, 

6 > + (1 - 9)u^2S + {l-9)\q- - 29] + 0{S^). (B.7) 



The curve £2 is nearly linear for ^ of order 1. It can be parametrized in the 
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rescaled coordinates [Ci,C2,'0] (see (B.2)) by the map 

L2:M^M':6^[6,??(6),6] 



(B.8) 



where the second component (^2) of -^2 is the locally unique solution ^2 of the 
equation 

/^V(6,6,ei,g,5) = o (B.9) 

depending on ^1, q and S (the grazing is quadratic, hence h'_^ is monotone in 
^1). More precisely, the solution ^2 of (B.9) can be represented as a graph 77(^1) 
(dropping the dependence of on g and S) , and this graph can be expanded in 
S: 

viCi) = -^-^S + 5'r{Ci,q,SMi,q]. (B.IO) 



The remainder term of order 5^ in (B.IO) is at most linear in ^1 and q because 
77(0) = for {1 = g = and all S (the grazing point ^1 — ^2 — i/^ — q — 
satisfies (B.9) for all 5). 



The map Hi We expand the trajectory {Ci{t):C2{t),'il>{t)) of £"+ through a 
point 

Po = (6,o,?7(Ci,o),Ci,o) e A 
in time (■0(^) = Ci(^) because po £ 



ei,o 



+ 







t + 



The coefficient in front of is non-zero in its first component because we have 
quadratic grazing. The first component of the coefiicient in front of t vanishes 
by choice of £2- Thus, for ^1^0 < and |^i,o| of order 1 or less, the trajectory 
intersects the switching manifold {^i = 0} at 



UJ 



(l + 0{S')) + 0{S')C^,o + 0{S% 



(B.ll) 



The map Hi is defined by following backward in time. Thus, wc choose the 
negative sign in (B.ll). This gives rise to an intersection of the trajectory with 
the switching surface at a point (0, ^2,1,0) with the coordinate 



6,1 = (1 + 0{6')) - + 0{6') 

' UJ 



(B.12) 



The first term of the sum in (B.12) contains all square-root terms of ^i^o- 
Consequently the map Hi maps 



Hi :po = [6,o,^(Ci,o),Ci,o]e>C2 = [0,6,i,0]e>C3 
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-3 -2 -1 

e 

Fig. B.2. Implicitly given graph of so{6). 

where ^2,1 is given by (B.12). The dominant term in the expansion of pi with 
respect to the coordinate ^i_o of po is a square root. 



The map 112 For ^10 of order 1 (and, thus, ^2,1 of order 1) the flow E_ points 
away from the switching surface Hg in pi because condition (B.7) is not satisfied 
in pi (because 9 < 0). The initial point for the map 112 is Pi = (0, ^2,1, 0). The 
trajectory (^i(s), ^2('5), '?/'(s)) following E_ through this point satisfies 



Us) 
Us) 



s 

6,1 



(B.13) 



s + exp(-s) - 1 + u;6,i[exp(-s) - 1]5 + 0(5^ 



where s = t/6 is the rescaled time. There is a unique time s*(^2,i) > such 
that the trajectory hits the switching surface {6*^1 + (1 — 9)4' = 0} again. The 
expansion of s=i,(^2,i) in ^ is 



s*{C. 



2,ly 



so{0) + S,;^^^.,, + O{5^) 



9 + so{9) 

and so{9) is the unique positive solution of the equation 

(l-0)exp(-so)-l + so + ^ = O. 



(B.14) 



(B.15) 



Figure B.2 shows the graph of so{9), which can be determined implicitly by 
expressing ^ as a function of Sq. Since we assume that \9\ is of order 1 we 
cannot replace (B.14) by its expansion in 9. The denominator in the coefficient 
for 6 in (B.14) is positive because the solution sq of (B.15) is greater than —9 
for all 6' < 0. Consequently, 112 maps 



^2-Pl = (0,6,1,0) ^P2 



s*(6,i),6 



2,l5 



-9sJ^ 
1-9 



2,1 
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For small S, and |^i^2| of order 1 or less, the sliding condition (B.6) is not 
satisfied for points in the image of 112 because Sq{9) > is of order 1. This 
implies that the flow will switch to E_^. The projection P3 of the image along 
the stable fibres of 



P3 ■ P2 



2,1) 



i-e 



maps into the slow invariant manifold of £'+ and away from the switching 
manifold because < —0/{l — 9) < 1 and 5=^(^2,1) > 0. Thus, the point p3 is 
also in the half space ^1+ = {(C, -0) : ^6 + (1 - > 0}. 



The map U3 The map U3 maps the point of the slow invariant manifold 
Ai^ back to the Poincarc section (more precisely its intersection with A1+), 
£2, by following E^. This can be obtained in the same manner as the map 
IIi: denote the components of ps by [G,3) '^2,3, "^^s] (where ips = ^i^s because 
Ps e and denote ^-component of the trajectory E^_^p3 by (Ci(^); C2(^))- 

The time ^3 when the trajectory hits £2 has to satisfy the condition defining 

>C2 

r?(6(^3))=6(^3), (B.16) 

which defines ts and the point p4 = (Ci,4, 6,4, ^/'4) = (6(^3), 6(^3), 6(^3)) 
implicitly. We expand ^(t) and S(^) with respect to time f in to second 
order in t using (B.3) (all third-order terms have a coefficient of order 5^) and 
(B.IO), and insert the expansion of ^3 in S, 

t3 = h,o + 6h,i + 0{6'), (B.17) 

into (B.16) where ^3,0 and ^3,1 are the unknowns at zero and first order level. 
Solving (B.16) to first order in S gives 

^3 = ^ + ^-±^5 + 0{S') and, thus, 
6,4= 6,3 + — <^ + C'(52). (B.18) 

Consequently, the map 113 maps 

Ha : P3 = [6,3, 6,3, 6,3] ^ P4 = [6,4, r?(6,4), 6,4] e £2 (B.19) 
where 6,4 is given (B.18). 
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Fig. B.3. The Poincare discontinuity map for the case 6 > near the point pg (0 in 
the rescaled coordindate system (B.2)). 

The PDM Combining the expressions (B.12), (B.14) and (B.18) we obtain 
the composition Pbm = lis o Pg o 112 o Hi: 

^'DM({,) = r-,„. . nil (B.20) 

We have expressed the map Pdm as a map for the coordinate of the point 
on the curve £2- Given the coordinate ^1, the other coordinates of a point of 
£2 are ^2 = ''7(^1) and ip = where the expansion of 77 is given in (B.IO). The 
expression for Sq{9) is given impUcitly by (B.15). The remainder term 0(5^) for 
^1 < may contain small corrections to the constant and square-root terms in 
^1. However, the form (B.20) guarantees that for sufficiently small 5 (and, thus, 
singular perturbation parameter e) the overall return map Pg o Pdm changes 
from a piecewise linear map for the reduced model {e = 0) to a discontinuous 
map for the full model (e > 0). The size of the jump at the discontinuity is of 
order e in the original coordinates. Moreover, the slope of the map next to the 
discontinuity is infinity from one side. 



B.2 The case ^ > 



For 6* > the assumptions behind the validity for the asymptotic form of the 
PDM near a grazing-sliding orbit are satisfied for the singularly perturbed 
system [4]. The flow Ej^ is grazing the switching surface Hs = {6x1 + {l — 6)y = 
0} quadratically in the point Pg = [—1,0,-1]^. The other flow E_ points 
toward the switching surface near the grazing point Pg of E^. Thus, according 
to Theorem 8.1 of [4] we have a generic grazing-sliding bifurcation. In general. 
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the zero-time discontinuity map Pz (ZDM) is piecewise asymptotically linear, 
having the form (ignoring higher-order terms in v) 



Pz{Pg + v) ^Pg + 



V — 



F- 



F+ 

dhF- 



if dhv > 0, 
dhv ii dhv <Q 



(B.21) 



where e is small, Pg is the grazing point, dh is the normal of the switching 
surface, and F± = {f±,g/e) is the derivative of E± in pg. In our example we 
have Pg = [-1,0,-1]^, dh = [^,0,1 -^], F+ = [0,-uj,Of and F_ = [1,0,0]^ 
at parameter = 1. Thus, at = 1 the linear approximation of the ZDM, Pz, 
for points in the invariant manifold of E+, is 



Xi 
Xi 



Xi 
X2 

■(i-r^)(xi + i)-i 

X2-uj9-\xi + 1) 

Xi 



if Xi > —1, 



if xi < -1. 



(B.22) 



The map Pz can be converted to the projected PDM Pdm by applying Pz to 
points in £2, then applying the projection P3 along the stable fibres of and 
finally applying the map 114, which is the same as for the case ^ < 0, following 
on its invariant manifold back to £2- Thus, ignoring terms of order 

3 /2 

Xi or higher, the PDM maps the coordinate xi of a point in £2 as follows: 



PDM(a;i) 



xi if xi > —1 

[1-^-^] {xi + l)-l iixi < -1. 



(B.23) 



The slope is different for xi < —1 from the slope of the reduced model {e — 0, 
where PDM(a;i) = —1 for xi < —1) if ^ e (0, 1). The difference in the slope 
is uniform for all small e > (expression (B.23) docs not depend on e). For 
parameters /i close to 1 the slope of Pqm for xi < —1 changes at most by a 
term of order ji — l. The slope of Pdm for xi < —1 is negative for 9 e (0, 1) 
and can be arbitrarily large in modulus, thus, giving rise to a strong expansion 
for small 9 (and, possibly, chaos) in a small neighborhood of the grazing orbit. 

The remainder of this section will show that the size of this neighborhood where 
the slope of Pdm is negative and non-small is at least of order (judging from 
Figure 4(b) it is larger). In order to see this we have to obtain the expressions for 
the maps Hi, 112 and lis for this case as well. The maps Hi and lis are identical 
to the case < apart from a slightly differently scaled set of coordinates. The 
map 112 follows the sliding flow Eg to the curve {/i'_,_(a;, y) = 0, h{x, y) = 0} in 
the case ^ > 0, and is, thus, different from the map 112 for the case ^ < 0. 
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We zoom into the neighborhood of the grazing point by introducing the new 
quantities ^, -0 and q with a scahng (different from the scahng for ^ < 0): 



xi = -1 + £ ^1, 
X2 = £6, 



// = 1 + e^q. 



(B.24) 



Similarly as for ^ < 0, any point in the phase-parameter space with ^, ijj and q 
of order 1 is in the vicinity of the grazing bifurcation. The switching manifold 
1-Ls is • ^^i + ll"^)"^ = 0}. System (13), (14) in the rescaled coordinates 

has the form (with rescaled time inew — toid/^ and truncating terms of order 
£3): 



-^6 + [el-26-?]£ + o(£3) 

ye 





if ^6 



if ^6 



{l-9)ij > 

{i-e)i/j < 



= -V' + 6 - ^6 + [^2 - 26 - q]e + 0{e') 



(B.25) 



(B.26) 



For system (B.25), (B.26) each trajectory can spend at most a time span of 
order 1 in any bounded ball (of size of order 1) in the (^ ^)-coordinates. The 
expansions of the quantities h'_^ and h'_ (determining where the system is 
sliding) in e are 



6 - + ki 



^.,-q-2^E + 0{6') 
e-'9 + ^i + {9-l)co^2 + 0{e). 



(B.27) 



As the expression for h'_ is always positive for ^ -0 and q of order 1 the system 
slides on Us whenever h'_^ < 0. 

The expressions for the curve £2 and the maps Hi and II3 arc identical to 
(B.8), (B.IO), (B.12) and (B.18), (B.19). The only modification is the differently 
scaled set of coordinates (B.24). The effect of this scaling is that we have to 
replace 5 by £ in all equations from (B.9) to (B.12) and in (B.18), leaving all 
expressions identical otherwise. Thus, Hi and map 

Hi : Po = [6,o,^(6,o),Ci,o] eC2^Pi = [0,6,1,0] €£3 
na : P3 = [6,3,6,3,6,3] ^P4 = [6,4,r/(6,4),6,4] e ^2 

where 77, 6,1 and 6,4 are given by 

26 + g 



^(6) 
6,1 



£ + £V(6,g,£)[6,g], 
26.0 + q 



(1 + 0{s')) - ^^i^£ + 0(s% 

UJ 



6,4= 6,3+^£ + 0(£2). 
UJ 
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The map 112 For ^i^o (and, thus, ,^2,1) of order 1 the flow points toward 
the switching surface "H^. In pi the flow points toward the switching surface, 
too (as /?/_!_ < in pi), such that 112 is defined by foUowing the shding fiow £"* 
until h'j^ (given by (B.27)) becomes zero. The shding flow El is invariant on the 
switching manifold Tig- It follows an ODE for ^ {ip is given hj t/j — —9 /{I — 6)^1). 
The right-hand-side of the ODE can be expanded in £ as 

e,= 4[&-&(i-»)l + o(.) ^^^^ 

6= -u + 0{e). 

The initial condition is (0, ,^2.1) (where ,^2.1 > 0). The time t that the system 
follows (B.28) until h'j^ becomes zero is of order 1 such that the trajectory ^[t) 
is a perturbation of order e of 

Ut) = u;(l - 9) ([1 - eM-t/0)] {Ouj + 6,1) - ^t) .3 29) 

We insert the trajectory (B.29) into the expansion (B.27) for h'^ to find the 
unique positive time ^2 for which e~'^h'_^ in (B.27) becomes zero in a point 

P2 = ^6,2,6,2, Q ^1,2 j ■ 

This time t2, the time for which for the system follows the sliding flow (B.28), 
is given implicitly to zero order in e by the equation 

- '"^ ^ + (l-^)exp(-t2/^) ■ ^^-^^^ 

The right-hand-side in (B.30) is uniformly monotone increasing in t2, and it 
is zero for t2 = 0. Thus, (B.30) can be solved for t2- In combination with the 
zero-order approximations 

^''^ - ^ (1 - ^) ^+(l-^)exp(-i2/^) (^-^'^ 
6,2 = 6,2/cu (B.32) 

for ^12 and ^2,2 one obtains well deflned graphs for ^1,2(^2,1) and ^,2(^2,1)- As 
monotonicity persists under small perturbations these graphs remain well- 
deflned also if we take all orders of e into account. 

When we combine the graph ^i,2('C2,i) with the maps Hi, the projection P3 
and the map 113 we obtain an implicit expression for the graph of the Pdm for 
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Fig. C.l. Dry friction oscillator: (a) setup, (b) friction force Fc between belt and 
block depending on relative velocity v^ei = x — v = X2 — v, (c) bifurcation diagram 
as predicted by simple model (C.l) with a grazing-sliding bifurcation at GS. 



^1 < 0, parametrized by the sliding time t2 G 



6,0 

6,4 



[1-9) [l-expi-t2/e)]+ t2 

e + {i-e)exp{-t2/e) 
e-[e + 12] exp{-t2/e) 



1 2 



'1 



exp(-t2/^) 



(B.33) 
(B.34) 



Its derivative at ^i^o = (that is, ^2 = 0) is {6 — l)/0. This is in line with the 
piecewise linearization (B.23), and gives evidence that e"^ is really the smallest 
present scale. 



C Mechanics example — autonomous dry-friction oscillator 



We demonstrate the practical consequences of Lemma 2 for a simple mechanical 
system: the autonomous dry-friction oscillator as shown in Figure C.l (a), a 
block of mass M attached to a rigid wall via a spring of stiffness K and slipping 
(or sticking) on a belt moving with velocity v toward the wall. This setup is a 
caricature of the mechanics of, for example, violins or door squeal. A simple 
model for the dry-friction oscillator is 

Xi = X2 

K 1 (C.l) 

^2 = -^xi--F,(x2-.;) 

where Xi is the position of the block, X2 is its velocity and is the friction 
force on the block exerted by the belt and acting against the relative velocity. 
The switching between vector fields is hidden in the form of the friction force 
Fc, which is the sum of a smooth nonlinear function and a term — /osign/i 
where h is h{x) = v — X2 for (C.l). See Figure C.l(b) for a typical model of Fc. 
Figure C.l(c) shows the bifurcation diagram predicted by model (C.l) (and 
qualitatively confirmed in experiments in [11]), showing that the equilibrium 
X = {—Fc{—v)/K, 0) loses its stability in a subcritical Hopf bifurcation at H 
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— identity 
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(a) e = 0.1, vq = 3.32 



4.1 4.2 4.3 4.4 4.5 4.6 
(b) £ = 0.01, vq = 3.1 



Fig. C.2. Approximate return maps P for the half-plane X2 = 0, xi > for dry 
friction oscillator (C.3). Other (non-dimensional) parameters are M = i^ = l,7 = 2, 
Fc{i^) = /osigni/ fyiy + fc{l- e"!'"!) signi^ where /o = -2, = 0.3, and fc = 2.5. 

and that slip-stick oscillations exist and are stable up to the belt velocity at 
point GS. The slip-stick oscillations correspond to periodic orbits with sliding 
segment (physically corresponding to sticking: the block has zero velocity 
relative to the belt for a while). At GS the slip-stick oscillations lose their 
stability (and cease to exist) in a grazing-sliding bifurcation as introduced in 
Section 4 and discussed for the simple minimal example in Section 3.1 (except 
that the orbits lying entirely in 71+ are unstable). Model (C.l) can be affected 
by singular perturbation in various ways. One possiblity is that the velocity v 
is not perfectly attained as a direct control parameter but instead a controlled 
motor is aiming to achieve a certain belt velocity vq- This would imply that v 
is determined as a dependent variable of a fast subsystem, for example 



My,eltV = a[vo-v]+ Fc{x2 - v) 



(C.2) 



where a is the proportional control gain and Mbdt is the effective mass of the 
belt. If we introduce the variable w = (M/Mbeit)a;2 + v and the parameter 
7 = M/Mbeit, and treat Mbeit/o as the small parameter e then system (C.l), 
(C.2) is in the slow-fast form (3), (4) as described in the introduction. Namely 
we have 



Xi = X2 
X2 



(C.3) 



ew = -w + vo + 'yx2 - e—xi. 



The slow variable is x = (xi,X2), the fast variable is w, and the switching 
function h is h{x, w) = w — {1 + ■j)x2- 
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The slow manifold Aio {e = in the fast subsystem of (C.3)) is the graph of the 
function w{x) = Vq + 7x2}. The first-order approximation of Ai± is the graph 
of w±{x) = vq + 7x2 + eFc{x2 — Vq). Thus, at the switching surface we have 
that w±{x) = Wo + 7a; =F ^fo/M, which means that one effect of the imperfect 
control (the non-zero e) is a slow-down of the belt velocity of order e in Ti+. 
The other, less obvious, effect is that the non-zero e can limit the minimal 
length of sliding (sticking) segments that we can observe for the stable stick-slip 
oscillations. Figure C.2 shows the approximate return map to the half-line 
£ — {{x,w) : w — w^{x),X2 = 0, xi > 0} close to the grazing-sliding point 
GS shown in Figure C.l(c) (again, we have neglected the exponentially small 
difference between w and w+{x)). The family of sliding periodic orbits (the 
stick-slip oscillations) undergo a saddle-node bifurcation and become imstable 
before the sliding segment shrinks to zero length. This is a consequence of the 
change of the slope in the graph of the return map in Figure C.2, a change which 
is of order 1 uniformly for £ — > 0. Apparently, the grazing-sliding bifurcation in 
the dry-friction oscillator falls into case 2 of Lemma 2. Namely, we observe a 
grazing-sliding bifurcation and the periodic orbit persists, but the stability of 
the periodic orbit may change completely. The mechanism behind this change 
is identical to the case ^ > of the minimal example in Lemma 3. We note that 
the change of the right-side slopes of the return map in Figure C.2 depends 
on 7 such that the saddle-node bifurcation occurs only for sufficiently large 7 
(that is, if the belt is not much heavier than the block). 
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